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We construct and compare a variety of simple models for strange stars, namely, hypothet- 
ical self-bound objects made of a cold stable version of the quark-gluon plasma. Exact, 
quasi-exact and numerical models are examined to find the most economical description 
for these objects. A simple and successful parametrization of them is given in terms of 
the central density, and many differences among the models are explicitly shown and 
discussed. 
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1. Introduction 

We study in this paper a variety of models for strange stars, namely compact objects 
constituted of an hypothetic cold version of a quark plasma featuring (by hypothesis) 
a lower energy per baryon than normal hadronic matter. Quarks were suspected to 
constitute a substantial fraction of a "neutron" star interior since the early '70pH2 
and their main properties explored. In 198421 the possibility of a star fully com- 
posed of quarks up to its surface (but perhaps covered with a "normal" crust) was 
considered. The simplest model was constructed with a free quark+vacuum energy 
equation of state (see Ref. 2}. Because the scale of the strong interactions is set 
by a multiple of the parameter B of the MIT Bag Model, identified as the vacuum 
energy density and numerically valued around 60MeV jmT z or so^, these compact 
objects feature masses ~ 1M© and radii ~ 10km, pretty similar to their "normal" 
cousins. In fact they may constitute all the "neutron" stars hitherto observed. 

The strong interactions (not gravity) bind these kind of compact objects, and at 
least for the low-mass range M < 1 M Q , are the responsible for the striking behavior 
of their mass-radius plot, namely the continuity of models all the way down to pretty 
small masses ~ 0.5M Q , the minimum allowed by the MIT equation of state. Even 
though the low-mass stars are clearly newtonian, gravity is increasingly strong at 
the higher end and forces the use of general relativistic models when ^f|jr > 0.2. 
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As is well-known, the models of these objects can be described, in a first approx- 
imation, by the solutions of the Einstein equations for a perfect isotropic, static and 
spherically symmetric fluid (whose metric is given by ds 2 = <?e vi - r ^dt 2 — e x ^dr 2 — 
r 2 d6 2 -r 2 sin 2 {6)dcj) 2 ): 

8ttG , . \'e- x l-e- A ( r ) 

p( r ) = __ + } (i) 



8ttG , . v'e~ x 1 - e- A ( r ) 

— P(r) = — , (2) 

P'(r) + ^(c 2 p(r)+p(r)y(r)=0, (3) 

where the last equation (contracted Bianchi identity) expresses the conservation 
of energy-momentum of the fluid. 

Mathematically, there are at least three different strategies to solve these equa- 
tions 

• if the pressure, or the density or one metric element is given (for example, 
by an ansatz), an exact or numerical solution can be found by integration. 
However, this does not guarantee to have any control over the equation of 
state p(p), and similar solution functions 

• instead, if the equation of state p(p) is given (i.e. the fluid is characterized 
from the beginning), the integration can be performed (at least numerically) 
and the properties of the stellar models follow 

• if both the equation of state and one additional function (p, p or one of the 
metric) are known, a match of the overdetermined system can be attempted, 
being in general possible for some values of the involved parameters only. 

It is not difficult to envision that the third route can be also employed without 
an overdetermined system since one provide more degrees of freedom (for example, 
electric field or pressure anisotropy) at the expense of modifying the originally posed 
physical problem. We shall return to this point below. 

In this work, we have performed three sort of studies: 

• the modeling of a strange star by using some of the well known exact 
solutions (only that ones that meet physical requirements) for the Einstein 
Field Equation for static and spherically symmetric perfect fluid, e.g., the 
solutions obtained by the first strategy; 

• the modeling of a strange star by the third routine, e.g. overdetermining 
the system, in order to not only have the microphysics involved (by means 
of the equation of state), but also to have control over the density profile, 
chosen to be gaussian as a first reasonable guess; 
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• a comparison of the above studies to the true solution, or profiles, for a 
strange star. By true solution, we mean that one obtained closing the system 
with the equation of state of MIT Bag Model (see below) and numerically 
integrating the system. 

With this we show that the applicability of the well known analytical solutions 
are rather limited for modeling these stars and the best one can do having in mind 
the original system is to construct a analytical quasi-exact solution with a overde- 
termined system, by the expense to introduce an error in the geometry. 

Then we show how the introduction of more degrees of freedom (specially pres- 
sure anisotropy and the existence of an electric field) can avoid an overdetermincd 
system and at the same provide self-consistent analytical solutions for strange stars. 

2. Calculations 

2.1. General considerations 

As a benchmark for our calculations we have adopted the simplest numerical inte- 
gration of the massless non-interacting quarks presented, for example, by Alcock, 
Fahri and Olinto (hereafter AFO) (Ref. |6]). This is mainly motivated by the shapes 
of the profiles obtained by them, which are reproduced in figures [TJ [2] and [3] It is 
clear that any other model in which interactions or finite-quark masses are needed 
can be obtained analogously (see Ref. 7 for an exploration of the full parameter 
space of the models). 

2 5 1 1 1 1 1 
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Fig. 1. Mass-central density relation of the strange star sequence calculated by Alcock, Farhi and 
Olinto Ref. 6 

Since the main goal of this work is to simplify the description of the stellar models 
as much as possible, possibly neglecting small (< 1%) differences arising from the 
fact that numerical and exact (or quasi-exact) models do not match exactly, it is 
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Radius [km] 

Fig. 2. Mass-radius relation of the strange star sequence calculated by Alcock, Farhi and Olinto 
Ref. 6 
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Fig. 3. Density profiles of four selected strange stars calculated by Alcock, Farhi and Olinto Rcf. 
6 

useful to parametrize all the features in terms of a single input parameter, chosen 
to be the central density p c which must satisfy p c > AB/c 2 in order to be consistent 
with the description of the matter in the MIT bag model, where p = (c 2 p — AB)/3. 

To fix ideas, assume that f(x) = ag(x) + b is a solution of the density profile. 
If an ansatz is given for g(x), with that profile it is possible to obtain solutions for 
the pressure and the for the metric elements. Then we fit the parameters a — a(p c ) 
and b = b(p c ), so that the mass-radius relation of AFO is reproduced as accurately 
as possible. After this is achieved, we can proceed to derive all the features of any 
given model. 

However, not all exact (or even quasi-exact, see below) solutions describe a viable 
star: physically acceptable models should fulfil the following conditions, as discussed 
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in the work by Delgaty and Lake (Ref. UJ) : 

• e A (°) = 1 

To integrate the Einstein equations, two initial conditions are needed, gen- 
erally taken as m(r = 0) and p(r = 0), this is why this condition is so 
important to define the central region of the star. A small sphere of radius 
r = e has a length of 2ttc and proper radius e A / 2 e. A small circle around 
r = has a ratio of the length to its radius given by 2ne~ x / 2 . However, 
since the space-time is locally flat, this ratio should be for this small cir- 
cle around r = 0, equal to 2n. Thus, e A (°) = 1 and if r — > 0, then also 
m(r) -> 0, cf. Ref. M 

• — constant < 1 

A particle characterized by the geodesic with a given constant £ has an en- 
ergy e~ v l 2 £ relative to a locally inertial observer at rest in this space-time. 
Since e~ v l 2 = 1 far away from the object, £ is the energy a distant observer 
would measure if the particle was distant (i.e. the energy at infinity). Since 
g-w/2 > i f or a jj ther regions not at the infinit jf], the observer will measure 
a higher energy than £ in that regions. This extra energy is the energy a 
particle gains when it falls in the gravitational field of the object. Thus, 
e -v/2 > ! ^ e „(o)/2 < i ( c f. Ref. O 

• p(r) regular at the origin and positive definite 

A region around r = 0, small enough to be called "center of the star" 
and large enough to maintain the physical features of the fluid at constant 
(uniform) density, must have a positive value for the fluid to be real; 

• p(r) regular at the origin and positive definite 

The pressure at a central core which is approximated by a constant (uni- 
form) density must have a constant positive value, as a result of the prop- 
erties of real matter; 

dp — 

The speed of sound in the fluid should be smaller than the velocity of light. 

Ou r fir st approach to the problem is described below, using the well-known Tol- 
man I\ffiS and Buchdhal ffH which satisfy all properties 1-5 necessary for realistic 
models. 

2.2. Exact and quasi-exact solutions 
2.2.1. Tolman IV and Buchdahl I 

Tolman IV and Buchdahl I are two well-known exact solutions featuring mathe- 
matically simple expressions for static self-gravitating fluid spheres. Moreover, they 
are physically sound and fulfil tghe requirements listed in the former section. We 



a From the continuity condition at the boundary of the object, e"( r ^) = 1 — < 1. But 
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have tested the applicability of these simple models to the strange star problem, 
comparing them with the numerical results of AFO. A previous attempt can be 
found in Ref. [12] 

We chose these two solutions also because the density profile can be well fitted 
in order to reproduce the AFO density profile. 
Tolman IV 

This solution has been found by R. Tolman in his seminal paper (|10|) . To solve 
the Einstein equations, he made an ansatz e v v /2r — const, a method that rendered 
an exact integration of the problem. The functions thus obtained read 

A 1 + 2r 2 /A 2 



{l-r 2 /R 2 )(l + r 2 /A 2 )' 
e» = Q 2 {\ 



r 2 



- . 



A 2 J' 

8ttG 1 l + 3A 2 /R 2 + 3r 2 /R 2 2 1 - r 2 /R 2 



-p(r) 



c 2 PW A 2 l + 2r 2 /A 2 A 2 (l + 2r 2 /A 2 ) 2 

and 

87rG 1 1 - A 2 /R 2 - 3r 2 /R 2 

c 4 PIT) - A 2 l + 2r 2 /A 2 

From these expressions, an equation of state, a physical radius (boundary, r^) 
the constant Q and the total mass of the sphere are, respectively 



pip) = Pc - 5-f +5 — 



Vc r P Q \ c 2 c 2 ) 
~2 ' 0_ 2 + Pa , > 



with 



_R_/ A2\V2 
H ~ 3Va V R?) ' 

2 _ 1 - rg/tf» 
V l + 2r 2 b /A 2 ' 

M _c 2 r b( ^ (l-r 2 /R 2 )(l + r 2 /A 2 ) ^ 

M bound . cond . - 2Q l + 2r 2 /A 2 )' 



8vrG n R 2 + A 2 SttG R 2 - A 2 

—5TPC=3 n O AO > — J-Pc 



R 2 A 2 ' c 4 ^ A 2 i? 2 
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with A, R and Q three arbitrary constants. 

To match numerical strange star models as calculated by AFO with the Tolman 
IV solution, we seek to reproduce the mass and the radius of each model, that is, to 
determine A, R and Q so that the mass and the radius of the numerical calculations 
are reproduced. However, the central density could not be adjusted simultaneously 
(the later depends on A and R, which in turn are determined by the radius and 
mass). A relevant question is whether the resulting equation of state is compatible 
with the linear expression of the MIT Bag Model. We can easily check that the 
Tolman IV is not appropriate, first because it is not linear, but more importantly 
because that it depends on the A and R present in the metric elements. When 
A and R are adjusted to reproduce accurately the masses and radii of the AFO 
calculations, they happen to depend on the central density of the model. Therefore, 
the equation of state varies from model to model (from star to star), in spite that 
its functional form remains the same. This behavior is illustrated in Fig|4l 



7 
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Fig. 4. Variation of the equation of state with the constants A, Ft present in Tolman's models 
Ref. 10 



It is clear that an equation of state that varies from model to model is not useful 
to describe a compact star, since the equation of state should reflect the fundamen- 
tal properties of cold matter and should not behave in this way. 

Buchdahl I 

Much in the same way as in Tolman IV, we examined the Buchdahl JIH solution. 
The expressions are given by 
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A _ 2(1 + Q 2 ) 
2-CY 2 ' 



i v = A[{1 + Cr 2 f /2 + B^2 - Cr 2 (5 + 2Cr 2 )\ 



o n __2Cr_ _ 4CV(2~CV 2 ) _ 2-Cr 2 

o7TCr , , 2+2CV 2 (2+2Cr 2 ) 2 1 2+2Cr 2 



-p(r) 



8lrG 2 [3Cr yTTC^ - BCr £+l G : ] + 4BCrV2 - Cr*\ {2-Cr 2 ) 



[(1 + CV 2 ) 3 / 2 + BV2-Cr 2 (5 + 2CY 2 )](2 + 2CY 2 )r 

-, _ 2-Cr 2 
1 2+2CV 2 

where A, B and C are arbitrary constants. 

Two types of fits were possible: in the first one we fitted the mass and radius 
of each model, and as a consequence p c ^ pc AFO as before. We also attempted to 
fit both the mass and central density, leading to different radii, depending on the 
masses. However, we eventually came to the same problem of the non-linearity of 
the equation of state. Again the obtained equation of state depends on variable 
quantities and is not useful to model the quark matter. 

2.2.2. Quasi-exact gaussian model 

An inspection to the density profiles presented by AFO (Figure [3]) suggests a simple 
and popular parametrization of the density dependence. As a simple ansatz we as- 
sumed a gaussian profile, a procedure that led to an overdetermination of the system 
(since the linear MIT equation of state was also imposed). With those assumptions, 
the following expressions were immediately obtained in closed forms 

P(r) = ( Pc -§)e-^ + §, (4) 
p(r) = c 2 p(r) - (5) 



"(0 = 23+ a. (6) 
_ x{r) _ 4 7 rG'p c e- r2 / r °rg _ 2-K 3 / 2 Gp c r z erf{r/r ) _ 4nGBe- r2 / r °r 2 
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2^/ 2 GBrlerf{r/r ) _ 8nGr 2 B 
+ re* 3c 4 ' ( ' 

r 2 (8irGr 2 e- r2 / r ° Pc c 2 - 8nGr 2 e~ r2 / r " B - 24irGr 2 B + 3c 4 ) 

3c 4( r 2 +r 2) ■ W 

The problem posed in this way reduces to prove whether the two expressions 
for exp(A(r)), eqs. (7) and (9), are equivalent. Alternatively, even if the functions 
are different, they could be almost identical inside the star, and therefore a small 
error would be introduced by employing either one. This is why we speak about 
quasi-exact solutions for the linear equation of state to stellar problem. 

We used as benchmark the same four density profiles shown in AFO for definite 
values of the total mass M = 1.99 M Q , M = 1.95 M Q , M = 1.4 M© and M = 0.53 
M . 

The idea was to fit expressions for a and for r a (the scale of decay of the density 
profile) as a function of the central density p c which could reproduce the masses 
and radii. 

After that, given any p c , one may obtain any stellar model along the sequence 
because he can calculate the total mass, the radius, the metric elements and the 
density and pressure profiles. 

However, it proved impossible to fit all parameters simultaneously without mak- 
ing the system inconsistent. Letting the radius to vary (i.e. relaxing the condition 
for them to reproduce the AFO results), but keeping the masses and central densi- 
ties, then the stellar radius R, the parameter a and the scale r a could be obtained. 
As an important remark, we stress that fixing the mass means that the integral of 
the density function and the boundary condition derived from the metric v{r) must 
be the same, with their values the same as the ones found by AFO. 

The fit resulting from just the four profiles did not allow a good determination 
of a and r Q . Thus, the grid of models was extended to 20 profiles (i.e. 20 values of 
the central density) to cover wider mass and radii intervals. The values of r Q and a 
are shown in TableTT] 

Having these values for a and r , we proceeded to derive an accurate functional 
fit of the form, 

a(w) = -2.4 + 9.96e(~^ + 2.426e (_ o 3 59) + 1.92e ( "^ ) (9) 



r (w) = 5.56 + 22.46e ( -^ } + 10.35e ( ~^ + 169.46e ( "o 3 r? ) [km], (10) 

where w = ^rl 2 p c , I = 10 km, is a dimensionless quantity of the order unity 
that embodies physical constants (c, G etc.) and the lengthscale of the object (~ 
10km). 
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Pc[W g/cm J 


M /Mq [Ar (J\ 


r [km\ 


a 


21 


2.014042750 


7.631010387 


-1.780521357 


20.5 


c\ a -i ,1 n a i ^ An 

2.014801608 


7.714337813 


i p?/" 1 n n n n i n 

-1.762999185 


20 


2.015326779 


*T AA1 A y( r AA1 

7.801045927 


^ PT/( /lAOAO ,1/7 

-1.744982347 


19.6 


2.015436678 


7.872850429 


-1.730140677 


18 


2.013326957 


8.185473708 


-1.666867564 


16 


2.002284428 


o f a f* i An Ann 

8.646108988 


1 rTTTATr AAA 

-1.577075908 


1 A 
14 


1.9 10 1 9/010 


9./lz9z9t)19 


1 /I 71 /I OOAQ 
-1.4 ( 14zzUo 


12 


1.921914333 


9.936363686 


-1.343215145 


11 


1.877755062 


10.38285995 


-1.267268177 


10 


1.815572219 


10.90882098 


-1.180584905 


9 


1.726746594 


11.54174779 


-1.079793495 


8 


1.597234612 


12.32531382 


-0.959831395 


7.5 


1.510155564 


12.79466558 


-0.890275006 


7 


1.402444810 


13.33296263 


-0.81247546 


6.5 


1.267766569 


13.95933096 


-0.724422383 


6 


1.097507111 


14.70148514 


-0.62334087 


5.7 


0.973622557 


15.21887892 


-0.554849121 


5.2 


0.721604642 


16.24225101 


-0.423820856 


4.885 


0.529204878 


17.02459731 


-0.327483126 


4.6 


0.333488500 


17.85691977 


-0.228380811 



With these values for r and a, we may wonder about the accuracy of the metric 
element e xl - r \ We have seen that an ambiguity for this quantity exists. Its first 
expression was derived from the first Einstein equation (eq. [TJ and the second from 
the second one (eq. [2]). How different they really are? Actually it can be checked 
that the agreement is very good for the lower densities (Figure [5]), and worsens 
for the highest densities (Figure HJ) along the stellar sequence. This behavior is 
a consequence of the newtonian character of the gravitational field for the lower 
densities, already discussed in Section 1. 

In summary, we have shown that with just the central density value (quite 
analogously to any other calculation of stellar structure that begins by specifying 
this value), the values of a and r a can be calculated from the analytical fit. They 
determine in turn the profiles of density and pressure within the gaussian ansatz, and 
with them the physical radius of the model R. Finally, the mass can be calculated 
analytically by integrating the gaussian density profile and/or using the boundary 
condition, yielding the same values within < 10% because of the deviations of the 
fitted expressions from the exact numerical results. We believe that this quasi- 
analytical model (eqs. 4-8) can be useful for a (very) accurate evaluation of the 
stellar model structure in a variety of situations. 

It is worth to mention at this point that Cheng and Harkcff^ presented a very 
similar approach, starting with the relativistic conditions of thermodynamic equi- 
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Fig. 5. The difference between the two expressions of the metric element e~ A ( r ' inside the star 
for the model with p c = 5.2 X 10 14 g/cm 3 . The radius of the star is marked by the vertical line 
at the end of the curves. Up: the two expressions; Bottom: the relative error as a function of the 
radial coordinate. 

librium and the MIT bag equation of state. They found very accurate approximate 
mass and radius formula for strange stars, in the static case and, amazingly, in 
the rotating case. Their errors are less than 1% for the former case and 3% for 
maximally rotating stars. 

Considering how the general relativity changes the conditions of thermal equi- 
librium inside the compact star, they started with the equation for the chemical 
potential 

dfi dp 
M e + p ' 

Comparing the conditions at the center and at the boundary (surface) 



Equation 7 




(11) 
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The same as in Fig. 5 for the model with p c = 19.6 X 10 14 <;/em 3 . The differences are 
as pointed out in the text. 



^ = ^_ fl 



2GM\ 
c 2 R ) 



A c{ Pc ,By > 1 (i2) 

and with the MIT Bag equation of the state, they could integrate the equation 
11 to obtain 



(13) 



/fg / ml/4 

{e c -B)^ {€ - B) ■ 

Here, the function C is related to the time component of the metric tensor at the 
center of the star, g^l c = C(p c ,B), that is also related to the redshift of a photon 
emitted at the center of the star. 

Thus, defining the dimensionless parameter 77 = p c /B, they obtained the follow- 
ing expression to the mass-radius ratio: 
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]. (14) 



M __ 
~R ~ 2G 

From considerations about the mass continuity on the surface, they proposed 
the following mass-radius relation 



M = ^-a(r,)R s . (15) 

Here, C(rf) — const + J^. ai/rf (a, constants) and similarly to 0(77). 

Finally, they integrated the Einstein Field Equations numerically and fitted 
C(j]) and 0(77) in order to obtain expressions for the mass and radius in closed 
forms, where they used B = 56MeV/ fm 3 . 

Their formulae, however, become increasingly inaccurate when p c — > AB. 

In the following, we give our derivation for mass and radius, the analogous to 
their formulae (see eqs. 14 and 15 of Ref. 13), 




We stress two important points about these two expressions. First, they pro- 
vide the "right" value for the zero-pressure situation, although this feature should 
not be taken too strictly, because the expressions are only approximations. Sec- 
ond, although they appear quite generic, one should remember that we used 
B = 57.5MeV/ fm 3 to derive a(w) and r (w). A different value of B would force to 
recalculate the latter. 

The advantage of our approach is that we have the profile for p(r), p(r), v(r) 

and A(r) in closed form for each possible value of the central density p c . Besides 

that, our formula are in complete agreement with the formula derived by Cheng 
m 01 

and Harko^ 21 . 

We also point out that recently, Narain, Schaffher-Bielich and MishustidHl have 
made a numerical calculation integrating the dimensionless TOV equation with 
a linear equation of state for generic fermionic matter. The possibility of scaling 
the solutions of the TOV problem was already known and exploited, for example, 
in Witten's paper Using this property, Narain, Schaffner-Bielich and Mishustin 
found a very useful general scaling solution and discussed how to rescale these 
equations in order to find the solutions for arbitrary fermion masses and interacting 
strengths. In spite of this generality, we point out that our process of making the 
equations dimensionless is different from that work and, moreover, the choice of 
a gaussian form induces the presence of the length r which controls the spatial 
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decay of the density. As a result, it is not possible to compare easily our quasi-exact 
models with their results, which remain more general but require the knowledge of 
the dimensionless M — R curve. 



2.2.3. Polytropic models 

A recent interesting approach has been presented by Lai and XtP^, which models a 
quark star with a polytropic equation of state. This kind of equation of state is gen- 
erally stiffer than that conventional linear one, like the bag model. Unfortunately, 
exact solutions are difficult to obtain for the polytropic case. However, their numer- 
ical results are interesting because they deal with two models: with and without 
QCD vacuum energy (A), and show how they compare to linear bag-like equations 
of state. See their numerical results in Figure 2 from Ref. [15] 

As a general feature we may say that the work of Lai and Xu shows that a 
polytropic description of strange stars is possible and accurate, but it does not 
provide an easy and economical form for generating models, at least not easier than 
a full numerical work. 



2.2.4. Exact anisotropic model 

The exact anisotropic model of Sharma and Maharaja is another accurate approxi- 
mation for the stellar structure and needs, as above, just one parameter (the central 
density). As a physical motivation to consider this anisotropy, they quote Usov^ 
that suppose a strong electric field formed into a thin layer at the quark surface of 
a bare strange star. This is possible due the depletion of s-quarks in this region. 
The exact expressions of this models are 



p(r) = 
m(r) 



6(3 + ar 2 ) 
8ir{l + ar 2 ) 2 '' 

r 3 b 



Pr(r) 



2(1 + ar 2 )' 
c 2 6(3 + ar 2 ) 4B 
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A(r) = (1 + ar 2 )c 2 
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t c 2 b(3+ar 2 ) jB\( Gr 3 b , 4nGr a p r (r) ^ - 
_ \ 127r(l+ar 2 ) 2 + 3 A 2 (l+aH) c 2 > ,„„.. 

I 2 Gr 3 b \ r v 6 ' 

and are enough to obtain the mass-radius relation, the density profile, and the 
pressures (tangential and radial). Here, b = 87rp c /3 and a = l/r„. We have checked 
that the anisotropy is small (see Figure [7} in the interval of central densities. Once 
the free parameters r and K are fitted as functions of the central density p c the 
models are completely specified. We show the resulting mass-radius relation in Fig- 
ure m 




ji 1 1 1 i 1 
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Radius [km] 




Radius [km] 



Fig. 7. Radial and tangential pressures for two anisotropic stellar models. Up: p c = 4.8 X 
10 14 g/cm 3 ; Bottom: p c = 24 X 10 14 g/cm 3 . 



The similarity of the curve in Figure [5] and the numerical results by AFO is ap- 
parent. However, in spite that masses and radii are practically the same, anisotropic 
models are systematically denser at the center by 25% or so. For example, the max- 
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Fig. 8. Mass-radius curve for the anisotropic model of Sharma-Maharaj Ref. 16 

imum mass model along the sequence has a central density of ~ 24 x 10 14 g/cm 3 
whereas in the AFO numerical calculation the value is just ~ 20 x 10 14 g/cm 3 . The 
density profiles are shown in Figure^ It is also apparent the similarities with AFO. 




Radius [km] 

Fig. 9. Density profiles of the anisotropic models of Sharma-Maharaj Ref. 16 

Rudermarff^l argued that the anisotropies could be important for densities p > 
10 15 g/cm 3 . This is consistent with our results. It should be pointed out, however, 
that the anisotropy could affect some parameters like the maximum mass (the effect 
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is small in our approach) and the redshift, as discussed long ago by Bowers and 
Lianfl 

Only the analytical solution for the anisotropic star is not sufficient for a com- 
plete description of such an object. This is because it is important to explain where 
the anisotropy comes from in order to produce the differences we have shown. Mak 
and Harko^ pointed out that a source of anisotropy could be an anisotropic veloc- 
ity distribution of the particles inside the star due, for example, a magnetic field, 
turbulence or convection. Perhaps the biggest challenge of the anisotropic model is 
the stability criterion. Chan et aPE showed that, in the onset of instabilities, even 
small anisotropics might drastically change the stability of the system. 

We have shown that the pressure anisotropy is small for low central densities 
but becomes larger and larger as the central density increases. Further studies are 
necessary in our approach to verify if all the sequence (e.g. in all central density 
range) is stable. 

2.2.5. Exact electric field model 

In the same way as before, we explored the model developed by Komathiraj and 
Maharaja to model strange stars parametrized by a single parameter, the central 
density. This model also has all the desirable properties cited above. However, the 
electric field is an explicit function of the position coordinate, starting from zero at 
the center and growing up to the surface. The effect of this field on the mass-radius 
relation is to increase the masses and respective radii. The mass-radius relation is 
shown in Figure 1101 
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Fig. 10. Mass-radius curve for the electric model of Komathiraj-Maharaj Ref. 22 
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In this model, we deal with a charged strange star. It is important to notice 
that this is quite different from the assumption made by Usov, since the charge is 
distributed inside the whole star. Charged compact stars also have been studied in 
many ways by Ray et a/Pi In this work they have found that a star can have a 
electric field about ~ 10 21 V/m in a particular case of polytropic equation of state. 
It is interesting to notice that the Komathiraj-Maharaj exact solution with electric 
field also provides electric fields at boundary > 1.5 x 10 w V/m. 

3. Conclusions 

We have used a wide set of mathematical solutions (exact and approximate) in 
order to model strange quark stars. All these solutions were parametrized by a 
single quantity, the central density. Both simple exact and quasi-exact models were 
addressed, with and without extra degrees of freedom (anisotropy and electric field). 

No simple analytical model was found to describe these stars by keeping the 
equation of state functionally unchanged (see 2.2.1). Solutions other than the ones 
discussed here should be examined having this description in mind. 

In our particular treatment, the quasi-exact solution, we have used a grid with 
20 points, corresponding to 20 central densities and we have fitted expressions for 
the free parameters of the problem that solves analytically the Einstein equations 
for a perfect fluid with spherical symmetry and a linear equation of state. The suc- 
cess of this particular approach can be gauged in Figure II H that summarizes the 
mass-radius relations of all solutions found and also in Figure [T^] where we show 
the density profiles of the stars of maximum mass from some solutions (Tolman, 
Buchdahl, gaussian-approximated and anisotropic). We remind that this work used 
the simple numerical models of AFO as a benchmark, but the later has a small 
number of degrees of freedom, being based on the MIT bag equation of state. This 
linear behavior, however, makes the fluid equation system impossible to be inte- 
grated analytically unless we give an extra equation overdetermining the system. If 
this is done, then it is necessary a match of the overdetermined system. 

We have shown that it is possible to find a quasi-exact solution that is useful to 
model a strange star in an easy way. We just give the central density and all other 
quantities like de mass, radius, profiles etc. can be found. With these analytical 
expressions we can predict many properties of the star defined by that central 
density. The errors remain < 1% in all quantities, including the geometry (Fig. 5 
and 6). 

We can see in the Figure [IT] that the gaussian-approximated solution is valid in 
the range of low central densities, matching almost exactly to AFO. The anisotropic 
model is valid for all range of central densities and is not very different of AFO. 
This is a promising model (see Figure |9] and compare it with Figure [3] but see also 
the remarks in the section 2.2.4). The electric field model is also valid, but only in 
the context of charged quark stars. 

This integrability problem does not happen with a system of fluid equations 
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featuring anisotropic pressure or an electric field. These new degrees of freedom 
makes the system integrable. We have discussed how these results compare to the 
linear, isotropic, uncharged models and how they affect the actual stellar features 
in practice. 

It should be reminded that, in spite of their deep physical differences, it is too 
early to dismiss any equation of state for strange stars. The extraction of the ra- 
dius from the observation of thermal-like emission is still problematic and may 
contain substantial inaccuracies. Taken at face value, however, the determinations 
have rendered small numbers, certainly incompatible with neutron star models if 
true. For example, the objects EXO 0748-67^231 [2.10 ± O.28M , > 13.8 ± 1.8km] 
and EXO 1745-248^1 [1.4M Q , llfcm] or [1.7M ,9fcm] can be described by the 
anisotropic model (taking into account the error bars) respectively [2.0M©, 11.44/cm] 
and [1.39M©, ll.Ofcm] (this last one with a claimed high precision). We conclude 
this work pointing out that simple, economical descriptions of a variety of strange 
stars can be constructed and may be useful to explore their properties and structural 
behavior in many static and dynamical situations. 
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Fig. 11. Mass-radius relations for various models. 
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